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ABSTRACT 

Disk galaxies are in hydrostatic equilibrium along their vertical axis. The pressure 
allowing for this configuration consists of thermal, turbulent, magnetic and cosmic ray 
components. For the Milky Way(MW) the thermal pressure contributes ~ 10% of the 
total pressure near the plane, with this fraction dropping towards higher altitudes. Out 
of the rest, magnetic fields contribute 1/3 of the pressure to distances of ~ 3kpc 
above the disk plane. In this letter we attempt to extrapolate these local values to 
high redshift, rapidly accreting, rapidly star forming disk galaxies and study the effect 
of the extra pressure sources on the accretion of gas onto the galaxies. In particular, 
magnetic field tension may convert a smooth cold-flow accretion to clumpy, irregular star 
formation regions and rates. The infalling gas accumulates on the edge of the magnetic 
fields, supported by magnetic tension. When the mass of the infalling gas exceeds some 
threshold mass, its gravitational force cannot be balanced by magnetic tension anymore, 
and it falls toward the disk's plane, rapidly making stars. Simplified estimations of this 
threshold mass are consistent with clumpy star formation observed in SINS, UDF, 
GOODS and GEMS surveys. We discuss the shortcomings of pure hydrodynamic codes 
in simulating the accretion of cold flows into galaxies, and emphasize the need for 
magneto- hydrodynamic simulations. 

Subject headings: galaxies: high-redshift — galaxies: ISM — galaxies: magnetic fields 



1. Introduction 

This letter deals with the accretion of the cold flows into disk galaxies at z = 2 — 4. These 
star forming galaxies are expected to have large magnetic fleld, to a degree that will signiflcantly 
affect the interaction of the cold flows with the disk. In particular, magnetic tension will cause 
smooth cold flows to enter the disk in a clumpy manner. We derive analytically the typical masses 
of these clumps and flnd that for reasonable parameters of the magnetic fields and the disks, the 
clump masses are 10'' — 10^ Mq as observed. The introduction also contains observational results 
regarding magnetic fields in the Milky Way (MW) and nearby galaxies (§ 1.1), and clump clusters 
and clumpy disk galaxies at redshift z = 2 — 4 (§ 1.2). In § 2 we derive a toy model for the magnetic 
tension of a disk (§ 2.1), typical parameters of filamentry cold flows (§ 2.2), and the predicted 
clump masses (§ 2.3). The importance of magnetic tension for gas accretion calculations and the 
conclusions are found in S 3. 
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1.1. The Pressure of the Milky Way and Nearby Galaxies 

The gaseous disk of the MW is highly multiphased and contains molecular gas, neutral gas and 
ionized gas at various temperatures and densities. Since these components are observed in many 
wavelengths, and by many instruments, it is difficult to combine the observation into one coherent 
picture of the interstellar medium (ISM). The vertical structure of the Galaxy is constructed of 
neutral gas and molecules dominating near the disk's plane, to neutral gas and ionized gas domi- 
nating at vertical heights larger than z ^ Ikpc, with a considerable ionized gas density extending 
out to ~ 3kpc. A compilation of the various components is available in Ferriere (2001). Using this 
information, one can compare the vertical gravitational attraction indicated by this column mass as 
well as stars and dark matter to the force arising from the thermal pressure of the combined gaseous 
components. It is striking that the thermal pressure at the mid-plane only supports ~ 10% of the 
pressure required to support the vertical extent of the disk (Badhwar & Stephens 1977; Dehnen 
& Binney 1998; Cox 2005). This fraction at mid-plane decreases further at larger radii. The re- 
maining 90 — 100% of the pressure at each altitude arises from a non-thermal pressure sources: 
turbulent, magnetic and cosmic ray pressures. The relative importance of the three components is 
comparable, indicating that they are roughly in equipartition (Boulares h Cox 1990; Beck et al. 
1996, and reference within.). This is currently a widely accepted hypothesis, which is motivated 
by the physical expectation that supernova explosions drives turbulence (McKcc & Ostriker 1977) 
which enhances seeded magnetic fields. The magnetic fields act to confine the cosmic rays and 
effectively set the diffusion rate of cosmic rays out of galaxies. 

Magnetic fields are measured by observing Zeeman splitting of 21cm, Faraday rotation of 
linearly polarised radio emitters and radio synchrotron emission of cosmic rays reacting to the 
magnetic fields. Radio loud nearby disk galaxies such as M51, M83 and NGC6946, have total field 
strengths of 20 — SOjuG in their spiral arms (Beck 2009, and reference within.) Bell (2003) showed 
that the radio emission, most of which is non-thermal synchrotron radiation, which arises from 
a rough equilibrium between cosmic rays and magnetic fields, correlates with SFR over 5 orders 
of magnitude in radio emission. The relation B ~ SFR°'^ is expected, and indeed observed, for 
star-forming galaxies (Lisenfeld, Voelk, & Xu 1996). We do not attempt to model this relation 
explicitly for our high-z star-forming (> lOOM© yr~^) galaxies, but adopt, as a lower bound, 
magnetic field amplitudes of 20 — lOO^uG as reasonable, conservative, extrapolation. A good review 
of the measurements and typical values of magnetic fields can be found in Beck (2009). 

The magnetic field of the MW, a low SFR, radio quiet galaxy, ranges from 6/iG at the solar 
neighborhood, to 10/iG in the inner Galaxy. At the spirals near the center, the magnetic fields can 
reach 100/iG. While the large scale structure of the MW's magnetic fields is hard to measure, the 
Galaxy provides a good probe of small scale magnetic fields. Magnetic fields seem to be entangled 
between the multiphase components of the ISM on parsec scales. 60% of the field's strength tends 
to be in coherent, ordered magnetic structures, correlating on kpc scales with the spiral structures 
of the MW. At higher altitudes the magnetic field drops, and at 3 kpc it is a factor of 2 lower when 
fitting synchrotron radiation, and five times lower assuming equipartition (Cox 2005). 
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It is worthwhile to emphasize the physical situation that arises from these conditions: the 
magnetic fields are pressurized but do not interact gravitationally. The field would disperse to 
infinity unless there is a small amount of gas locked between the entangled magnetic field, pulling 
it towards the disk. While the gas does not contribute a significant amount of thermal pressure, it 
contributes all the gravitational force, pulling everything together. 

1.2. Clumpy Star Formation at High Redshift 

During the past few years data have begun to accumulate regarding the structure of high 
redshift spiral galaxies. Their morphology was investigated using the Hubble Deep Field survey 
Cowie, Hu, & Songaila (1995) and the Ultra Deep Field (Elmegreen et al. 2005, 2007). They 
showed that the majority of high-z galaxies are disks, with clumpy, irregular morphology ("Clump 
clusters" or "Chain galaxies" ) and high velocity dispersions. Most of the star formation is observed 
to be concentrated in a few, well defined clumps with spatial extent of ~ Ikpc and mass of 10'' — 
10^ Mq. The inferred velocity dispersion of these galaxies is 20 — 50km sec~^. These observation are 
complemented spectroscopically by observations by the SINS survey (Genzel et al. 2006; Forster 
Schreiber et al. 2006) that contains 80 Ha emitting galaxies with masses M ~ 3 X 10^°Mq, most 
exhibiting rotating irregular gas patterns (Shapiro et al. 2008). Lower redshift information, derived 
from the GOODS and GEMS surveys, is analyzed in Elmegreen et al. (2009) and from the MASSIV 
survey in Epinat et al. (2009). Elmegreen et al. (2009) demonstrated that there is an evolutionary 
sequence of galaxies with the highly irregular clumpy ones at one end (making 50% of the galaxies 
at z ~ 4), and a smoother disk at the other (making 90% of the galaxies at z 1). 

While these galaxies are star-bursting, with typical SFR of 100 — 2OOM0 yr~^, most show no 
signatures of galaxy mergers and appear to maintain their disk. Their gas supply is most likely a 
result of gaseous cold flows accreting efficiently onto the galaxies (Dekel et al. 2009b; Bournaud & 
Elmegreen 2009) . Gas that is accreted onto the halo is unable to sustain a virial shock because the 
"would-be" post-shock gas would be cooling faster than it is contracting - unable to provide pressure 
support for the virial shock. This has been shown analytically by Birnboim & Dekel (2003); Dekel &: 
Birnboim (2006) and numerically by Keres et al. (2005); Ocvirk, Pichon, & Teyssier (2008); Keres 
et al. (2009); Brooks et al. (2009); Agertz, Teyssier, &: Moore (2009). There is a broad transition 
(Dekel &: Birnboim 2006; Keres et al. 2009) between cold-dominated halos and hot-dominated halos 
at which the accretion is in narrow, cold filaments penetrating through a hot, diffuse halo. The 
cold fiows (~ W^K gas) fiow undisturbed, preferentially perpendicular to the disk's plane (Dekel 
et al. 2009b), until the vicinity of the disk. These filaments then deposit their kinetic energy near 
the disfl join the disk and efficiently form stars. Based on lambda cold dark matter (ACDM) 



^The slowdown can be gradual, by a series of weak shocks (Keres et al. (2009); Agertz et al. (2009) using SPH 
and AMR simulations) or abrupt in a strong, "isothermal" shock (Dekel et al. (2009b); Ocvirk et al. (2008), AMR 
simulations) 
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hierarchical formation models Neistein & Dekel (2008); Dekel et al. (2009b) find that 2/3 of the 
accreted gas are in the form of smooth, cold, inflowing filaments. 

Recent theoretical work attempts to show how such clumps can be created. By using hy- 
drodynamic simulations of disks, Bournaud, Elmegreen, & Elmegreen (2007) demonstrated that 
gas rich disks corresponding to these high-z galaxies break rapidly into clumps via gravitational 
instability, and that these clumps migrate towards the center, creating a bulge. Dekel, Sari, &: 
Ceverino (2009a) showed that when there is a constant supply of smoothly accreted gas, the disk is 
maintained marginally Toomre stable, allowing for a prolonged evolution of bulge and clumpy disk 
consistent with observations. Using high resolution AMR simulations, Agertz et al. (2009) found 
that clumps naturally form within the disk from cold flows, and that the cold flows do not exhibit 
strong shocks at the edge of the disks. 



2. Disk Accretion in the presence of Magnetic Tension 

2.1. Magnetic Tension 

In this section we estimate the dynamic effect of magnetic tension on the accretion efficiency of 
cold-filament gas onto the galaxies. If one pulls a magnetic tube along its field lines, the amplitude 
of the magentic field does not change, while the volume it occupies grows. The work required to 
extend the tube by a distance dl is: ^Adl with A being the surface crossection of the tube. The 
tension of the magnetic field is thus: 

^ = ^. (1) 

Sir 

The magnetic fields are entangled on some scale, which creates a net of magnetic fields woven 
through the gas. While the field adjusts quickly to pressure variations, it does react against being 
bent. The toy model proposed here treats the field as a sheet of elastic matter, being bent inwards 
by gas that is accreted onto it. As gas falls on top of the magnetic sheet, it will fiow and settle at 
the lowest altitude it can, to minimize its potential energy. We assume that the deflection of the 
sheet is cylindrically symmetric around the center of the perturbation, and that the magnetic fields 
are pinned to their initial altitude at some radius by the joint forces of the entangled magnetic 
fields. We denote the radial direction by r and the vertical dimension by z. The magnetic fields 
are pinned to z^ax at rmax (fig- [I])- We solve for the shape of the line z{r) for which the system is 
static in the presence of mass Mc\ of gas with density Pgas- The gas will fill the well in the magnetic 
fields from the lowest z possible to some height Zgas according to: 

rr(Zgas) 

Mcl = / 27rr/)gas(2:gas - z{r))dr. (2) 

JO 

The force on some ring element of gas between [r : r + Ar] is the sum of the gravitation on that 
element, and the difference between the vertical component of the tension at r + Ar and r: 

2vrr(zgas - z{r))p^^g{z) dr = (3) 
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27r(r + Ar)Wasm[9{r + Ar)] - 27rrWasm[e{r)], 

with g{z) being the gravitational acceleration caused by the dislsEI, 9{r) the angle between the 
direction of the line and the horizontal direction, and W the width of the elastic layer (see fig. [1]). 
Throughout this work we assume that the width of the magnetic layer being bent scales like the 
amplitude of the perturbation, up to a factor 

W = lw^ - ^(0)). (4) 

A more sophisticated decision would require observational data or MHD calculations about the 
shape of the magnetic fields. The solution of eq. [2] and [3] is found by integrating from to rmax with 
[-Zgasi -2(0)] as boundary conditions, and fitting them to find the required values of [M^, z(rmax)]- 
by the simplex numerical method (Press et al. 1996). 



2.2. Magnetically Supported Gaseous Blobs 

In an effort to relate this magnetic tension to realistic cosmological conditions some typical 
parameters are chosen. The following is a rather lengthy "back of the envelope" calculation intended 
to estimate the density of the infalling gas as it interacts with the magnetic field. The puncture 
mass, actually, does not strongly depend on this density, but it is of deductive value to work out 
these numbers non-the-less. Based on simulations we assume 3 conical filaments, each at roughly 
20 = 20° in diameter. The temperature of the filaments is Tgi = 2 x IQ^K and they are in 
hydrostatic equilibrium with the virialized hot gas which is at T^ir- The density of the filaments in 
these conditions is p^i = fhPu^OD[iT^\/T^K){X ~ 37r0^/47r) + 30^/4]"-"^, with fh-,Pu and Sq-q being 
the universal baryonic fraction, density and virial overdensity respectively. We assume that the 
velocity of the filament is roughly the sound speed of the hot component (c^ = ^ksT^uNA/ p)- 
The global accretion rate M for halo mass M^i^ at redshift z is taken from Neistein &: Dekel 
(2008); Dekel et al. (2009b) which sets the radius of each of the three filaments near the disk 
according to: M/3 = -kB^yCsPsi- Finally, the infalling gas in the filament is stopped, converting its 
kinetic energy to thermal energy. The strong shock or set of weak shocks would adjust themselves 
so the ram-pressure of the infalling gas is balanced by the dense, halted gas below. Assuming 
that the post-shock pressure is all thermal the filament density at the edge of the disk is Pedge = 
Pfuc'^/ikBTfiiNA /i^^). For Mhaio = W^^Mq,z = 2 we have: p^dge = 8 x 10-24gr cm'^, = 



^throughout this work we use a fit for the vertical gravitational profile from Cox (2005) that takes into account the 
gas components, as well as the star and dark matter contribution. For different disk column densities the gravitational 
force is scaled, but its shape and z-dependence is not changed. The observed values of column densities by Elmegreen 
et al. (2009) of z ~ 2 - 4 galaxies (~ lOOMopc"^) are roughly the same as those of the MW. Since a and g in eq. 
[3] are proportional, a change in the scaling of the gravitation corresponds to a change in the square of the magnetic 
fields. 
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180km sec = TSM© yr ^, Rqi = 7kpc and for z = 4 we have: Pedge = 7 x 10 ^^gr cm ^, = 

230km sec-^M = 25OM0 yr"i and Rm = 4kpc. 

The final densities at the edge of the disks are upper limits. If some of the pressure of this 
gas is in turbulent form, or if the energy carried by the filament is spread over a larger volume, 
this density would be lower. It can be easily shown that the gas at the edge of the disk is in rough 
pressure equilibrium with the total pressure of the disk as is the case for the MW. However, while 
the MW's density is smooth and monotonically decreasing as z increases, no monotonic density 
profile is possible for cold accretion. When the temperature above the disk is similar to that of the 
disk, the extra pressure terms at the disk require that the density above the disk be larger than 
at the disk. This gas would sink in (via Rayleigh Taylor instabilities) unless the disk has surface 
tension such as the magnetic tension discussed here. The ionization fraction of the infalling gas is 
also of interest because the magnetic fields interact only with ionized gas. We expect that gas at 
1 — 2 X W^K and densities as mentioned above the gas will always be sufficiently ionized so that 
ambipolar diffusion will be ineffectiv^. 



2.3. Typical Masses of Magnetic Tension Induced Clumps 

The mass of clump for which the clump can puncture through the magnetic tension depends 
on the magnetic strength (eq. [3|), width of the sheared layer (eq. HI Ziy = 0.3. throughout this 
section) , the gravitational profile of the disk, the density (through eq. [2] and [3]) and on the radius 
at which magnetic fields are assumed to be pinned to their location by the global forces of the disk 
{Rmax = 3kpc throughout this section). Figure [2] shows four clump masses at density 10~^^gr cm~^ 
and magnetic field of lOO/xG. As the mass increases, the gravitation bends the magnetic fields and 
gas fills the lowest potential energy possible. The puncture mass is just above Md = W^Mq. 

Figure [3] shows the puncture mass for various gas densities and magnetic fields. For the MW, 
the magnetic tension will not affect the accretion of gas onto the galaxy (the "MW" rectangle 
in the figure). However, for larger magnetic fields and reasonable gas densities (marked with the 
"high-z" rectangle in fig. [3]) the mass of clumps that can penetrate is determined by the magnetic 
tension. For a wide range of parameters, particularly for high magnetic fields and lower accretion 
densities, the puncture mass increases to > 10^^ Mq. This gas will break through the magnetic fields 
instantaneously, yielding large gas fractions which are necessary for the in-situ clump instabilities 
found in Bournaud & Elmegreen (2009). The degree of coherence of ages of stars in the clumps is 
a good test for this scenario. 



■^The timescale for ambipolar diffusion of neutral gas with number density and ionization fraction of n = 
1 cm~^, X = W~'^ through magnetic fields B ~ 20/iG over a distance of ~ 1 kpc is t > 10^^ yr (Brandenburg 
& Zweibel 1994). 
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3. Discussion 

We have shown that for z = 2 — 4 star forming galaxies the magnetic fields at the edge of the 
disk are strong enough to alter the behavior of gas falling onto these disks. The expected magnetic 
tension is capable of stopping the cold accretion until some threshold mass is accumulated, at which 
point the chimp punctures through the magnetic fields and falls into the disk. The model could 
either produce individual clumps, or suppress accretion onto the galaxy for long periods of times, 
causing critical, coherent star-bursts on global galactic scales. Since the estimation of the magnetic 
field strengths here is very conservative, the magnetic tension will probably also affect gas accretion 
in more typical star forming galaxies with rates of < lOOMQyr"^. The model depends on the 
shape, correlation length and strength of the magnetic fields which are poorly known by theoretical 
considerations and extrapolation from local galaxies. There are no observational constraints on 
magnetic fields at these redshifts. 

For the MW, and nearby quiescent spirals, the magnetic fields are not strong enough to sig- 
nificantly affect the gas accretion, and the SFR will be dominated by the supply rate and other 
processes. Even if the magnetic fields were stronger, the time that would take for significant clumps 
to accumulate (at 1 — 4M0 yr~^) would be long with respect to the time it takes for magnetic fields 
to realign, which should be of the order of the dynamical time of the disk. This is not the case 
when the accretion rate is > IOOM0 yr~^. 

Neutral gas can potentially flow through the magnetic fields. The conditions required for such 
a process are typical for molecular clouds with extremely low ionization fractions. The relatively 
low densities and high temperatures (Tgas = 2 x lO^K) discussed here will render this ambipolar 
diffusion negligible. 

The toy model described in § 2 should not be taken at face value. It should, however, serve 
as an indication that non-thermal pressure components will have a significant, non-trivial dynamic 
effect on the accretion of gas onto galaxies. Pure hydrodynamic simulations are missing the effects 
of magnetic fields and cosmic rays, which account for ~ 2/3 of the pressure in the disk. This could 
be compensated in part by a correct choice of the cooling function or the equation of states of star 
forming regions (Springel &; Hernquist 2003). However, no form of gaseous equation of state can 
produce stress-strain behavior similar to that of magnetic tension. If this is important, full MHD 
codes should be used. Hydrodynamic simulations will predict correctly the shape and mass of the 
gaseous halo, the rate of the infalling gas and the cold filaments. They will also predict correctly 
the overall SFR of galaxies, because it is basically a question of gas supply. They are not expected 
to correctly predict the local, or temporal gas accretion or SFR, when the accretion is faster than 
typical mixing times in the disk. Other processes, such as angular momentum transfer in the disk, 
will depend on the dynamics within the disk which are dominated by non-thermal pressure. Shocks 
and sound waves in the disk will depend on the turbulence and Alfven speeds, rather than on the 
sound speed of the gas. 

I thank Donald Cox for spending time and effort to walk me through the magnetic field 
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calculations, and Shmulik Balberg for useful discussions and good advice. I thank the referee, 
Frederic Bournaud for useful comments. 
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Fig. 2. — The evolution of the magnetic layer as the mass is increased from 10^ — 10^ Mq. The 
magnetic fields react to the weight of the gas by bending inwards, until, around 10^ Mq the tension 
cannot support the weight any more and a static solution becomes impossible. The gas density, 
magnetic field, shear layer width and maximal radius of perturbation are pgas = 10~^^gr cm~^, B = 
100/xG, Iw = 0.3, -Rmax = 3kpc. 



- 12 - 



-21 
-22 
-23 
-24 
-25 
26 
-27 



f 107-' , 



/ 



/ 



/ 



/ 



/ 





10 



100 



Fig. 3. — The dependence of the puncture mass (various hues) on the magnetic field (x-axis) 
and density (y-axis). The rectangles mark typical values for the MW and high-z galaxies in that 
parameter space. 



